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1 Introduction 

Mixing transformations of fermion fields play a crucial role in high energy physics. 
The original Cabibbo mixing of d and s quarks and its extension to the Kobayashi- 
Maskawa three flavors mixing are essential ingredients in the Standard Model phe- 
nomenology [1]. On the other hand, although clear experimental evidence is still 
missing, it is widely believed that neutrino mixing transformations are the basic tool 
for further understanding of neutrino phenomenology as well as of solar physics [2]. 

In contrast with the large body of successful modelling and phenomenological 
computations (especially for Cabibbo-Kobayashi-Maskawa (CKM) quark mixing, the 
neutrino mixing still waiting for a conclusive experimental evidence), the quantum 
field theoretical analysis of the mixing transformations has not been pushed much 
deeply, as far as we know. The purpose of the present paper is indeed the study of 
the quantum field theory (QFT) framework of the fermion mixing transformations, 
thus focusing our attention more on the theoretical structure of fermion mixing than 
on its phenomenological features. 

As we will see, our study is far from being purely academic since, by clarifying 
the theoretical framework of fermion mixing, we will obtain some results which are 
also interesting to phenomenology and therefore to the real life of experiments. 

In particular, to be definite, we will focus our attention on neutrino mixing trans- 
formations and our analysis will lead to some modifications of the neutrino oscillation 
formulas, which in fact may be subject to experimental test. 

The paper is organized as follows. In Section 2 we study the generator of the 
Pontecorvo neutrino mixing transformations (two flavors mixing for Dirac flelds). 
We show that in the Lehmann-Symanzik-Zimmermann (LSZ) formalism of quantum 
fleld theory [3,6,7] the Fock space of the flavor states is unitarily inequivalent to the 
Fock space of the mass eigenstates in the inflnite volume limit. The flavor states are 
obtained as condensate of massive neutrino pairs and exhibit the structure of SU (2) 
coherent states [4]. In Section 3 we exhibit the condensation density as a function 
of the mixing angle, of the momentum and of the neutrino masses. In Section 4 we 
derive the neutrino flavor oscillations whose amplitude turns out to be momentum 
and mass dependent. This is a novel feature with respect to conventional analysis 
and may be subject to experimental test. In some sense, from the point of view of 
phenomenology, this is the most interesting result. Nevertheless, the condensate and 
coherent state structure of the vacuum is by itself a novel and theoretically interesting 
feature emerging from our analysis. In Section 5 we extend our considerations to three 
flavors mixing and show how the transformation matrix is obtained in terms of the 
QFT generators introduced in Section 2. We also obtain the three flavors oscillation 
formula. Finally, Section 6 is devoted to the conclusions. Although the group the- 
oretical analysis is conceptually simple, specific computations are sometime lengthy 
and, for the reader convenience, we confine mathematical details to the Appendices. 
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2 The vacuum structure for fermion mixing 

For definitiveness we consider the Pontccorvo mixing relations [5], although the 
following discussion applies to any Dirac fields. 
The mixing relations are: 

^'e(x) — i/ilx) cos^ + f2(x) sin^ 

u ij,{x) — —i'i(x) sin9 + U2{x) cos9 , (2.1) 

where I'dx) and ^'/^(a;) are the (Dirac) neutrino fields with definite fiavors. i'i{x) and 
1^2 (x) are the (free) neutrino fields with definite masses mi and m2, respectively. The 
fields (x) and 1^2 (a^) are written as 

= :h T.KA,i + ^"^'•']' ^ = 1' 2 ■ (2-2) 

^ k,r 

In the following, for simplicity, we will omit the vector notation for k and use the same 
symbol k to denote both k and its modulus k. al ■ and (31 ,-,2 = 1,2, r = 1, 2 are the 
annihilator operators for the vacuum state |0)i,2 = |0)i®|0)2: ttfc,i|0)i2 = Pk,i\^)i2 — 0- 
In eq.(2.2) we have included the time dependence in the wave functions. In the 
following this dependence will be omitted when no misunderstanding arises. The 
anticommutation relations are: 

{uri^),^Piy)}t=t' = S'{x - y)5^p5,, , = 1, ..,4 , (2.3) 

and 

Wk,i^ "gj} = hq^rs^ij] {(^li, (3fj} = SkqSrsSij, t,J = 1,2 . (2.4) 

All other anticommutators are zero. The orthonormality and completeness relations 
are: 

E ra* sa _ ST^ ra* sa _ r \^ ^.ra* sa _ \^ ra* sa _ q 

EK7< + vLrA) = ^"/^ ■ (2-5) 

r 

Eqs.(2.1) (or the ones obtained by inverting them) relate the respective hamiltonians 
Hi^2 (we consider only the mass terms) and iifg,// [5] : 

Hl^2 = mi + 7712 (2.6) 

He,u, = mee vl^e + "^^^ V^V^ + {ylv^i + i^Ji^e) (2.7) 

where = m\ cos^ 6 + m2 sin^ 6, 171^^ = mi sin^ 9 + m2 cos^ 9 and mg^ = (m2 — 
mi) sin ^ cos ^. 
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In QFT the basic dynamics, i.e. the Lagrangian and the resulting field equations, 
is given in terms of Heisenberg (or interacting) fields. The physical observables are 
expressed in terms of asymptotic in- (or out-) fields, also called physical or free fields. 
In the LSZ formalism of QFT [3,6,7], the free fields, say for definitiveness the in- 
fields, are obtained by the weak limit of the Heisenberg fields for time t —oo. The 
meaning of the weak limit is that the realization of the basic dynamics in terms of the 
in-fields is not unique so that the limit for t ^ —oo (or t ^ +00 for the out-fields) 
is representation dependent. Typical examples are the ones of spontaneously broken 
symmetry theories, where the same set of Heisenberg field equations describes the 
normal (symmetric) phase as well as the symmetry broken phase. The representation 
dependence of the asymptotic limit arises from the existence in QFT of infinitely 
many unitarily non-equivalent representations of the canonical (anti-) commutation 
relations [6,7]. Of course, since observables are described in terms of asymptotic fields, 
unitarily inequivalent representations describe different, i.e. physically inequivalent, 
phases. It is therefore of crucial importance, in order to get physically meaningful 
results, to investigate with much care the mapping among Heisenberg or interacting 
fields and free fields. Such a mapping is usually called the Haag expansion or the 
dynamical map [6,7]. Only in a very rude and naive approximation we may assume 
that interacting fields and free fields share the same vacuum state and the same Fock 
space representation. 

We stress that the above remarks apply to QFT, namely to systems with infinite 
number of degrees of freedom. In quantum mechanics, where finite volume systems 
are considered, the von Neumann theorem ensures that the representations of the 
canonical commutation relations are each other unitary equivalent and no problem 
arises with uniqueness of the asymptotic limit. In QFT, however, the von Neumann 
theorem does not hold and much more careful attention is required when considering 
any mapping among interacting and free fields [6,7]. 

With this warnings, mixing relations such as the relations (2.1) deserve a careful 
analysis. It is in fact our purpose to investigate the structure of the Fock spaces 7Yi,2 
and ?ie,^t relative to i^i(x), t'2(x) and I'dx), i^fiix), respectively. In particular we want 
to study the relation among these spaces in the infinite volume limit. We expect that 
7ii^2 and Ti-e^n become orthogonal in such a limit, since they represent the Hilbert 
spaces for free and interacting fields, respectively [6,7]. In the following, as usual, 
we will perform all computations at finite volume V and only at the end we will put 
V ^00. 

Our first step is the study of the generator of eqs.(2.1) and of the underlying group 
theoretical structure. 

Eqs.(2.1) can be put in the form: 

u^{x) = G-\e) <(x) G{e) 

M = G-\0) v^{x) G{e) , (2.8) 
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where G{9) is given by 



G{9) = exp 



(2.9) 



and is (at finite volume) an unitary operator: G~^{9) — G{—9) — G\9). We indeed 
observe that, from cqs.(2.8), cPv"/d9'^ — —v" , cPi''^/d9'^ = — z/^. By using the 
initial conditions 1^^10=0 = ^i, du"^ /d9\e=Q — 1^2 ^'^\o=o — ^2-1 di''^/d9\e=o — —1^1, 
we see that G{9) generates eqs.(2.1). 
By introducing the operators 

S+ = j d^xvl{x)v2{x) , S- = j d^x vl{x)vx{x) =(-5+)^ , (2.10) 

G{9) can be written as 



G{9) = exp[9{S+ - S-)] . 
It is easy to verify that, introducing 5*3 and the total charge Sq as follows 



1 f 

'S'3 = - y d^x (ul{x)iyi{x) - z/|(a;)z/2(a;)) , 



(2.11) 



(2.12) 



5-0=^1 d^x {ul{x)u,{x) + 4{x)u2{x)) , (2.13) 
the su{2) algebra is closed: 

[5+,5_]=253 , [53,^±]=±^± , [5o,^3] = [^0,^±]=0. (2.14) 

Using eq.(2.2) we can expand 5*+, S-, S3 and Sq as follows: 

k k r,s 

{'^k,l'^k,2 <^k,l^k,2 + '^-k,l'^k,2 P-k,l^k,2 + '^fe|l'^-ifc,2 <^k,lP-k,2 + '^-k,l'^-k,2 P-k,lP-k,2) ) 

(2.15) 



('^ifc!2'^M ^k,2'^k,l + '^-k,2'^k,l /^-fc,2'^fc,l + '^k,2'^-k,l (^k,2P-k,l + '^-k,2'^-k,l P-k,2P-k,l) ' 

(2.16) 



rt 



^4 



S3^^S', = -Y: - f3%f3L,,^ - ai«^,2 + P-i2P-k,2) , (2-17) 



2 k,r 
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So^Y: S'o = 2 E - P-k,iP-k,i + «^!2«^,2 - P-UP'-k,2) ■ (2.18) 

k k,r 

It is interesting to observe that the operatorial structure of eqs.(2.15) and (2.16) is the 
one of the rotation generator and of the Bogohubov generator. These structures will 
be exploited in the following (cf. Section 3 and Appendix D). Using these expansions 
it is easy to show that the following relations hold: 

[SlS'l]^2Sl , [Sl,Si]^±Si , [S^,Sl]^[S^,Si]^0, (2.19) 



[Si,Si] = [Sl,Si] = [S'„Sl]=0 , k^p. (2.20) 

This means that the original su{2) algebra given in eqs.(2.14) splits into k disjoint 
suk{2) algebras, given by eqs.(2.19), i.e. we have the group structure (Sik SUk{2). 

To establish the relation between Tii 2 and He,n we consider the generic matrix 
element 1,2(0!^" (3^)1 1,2 (a similar argument holds for z/^(a;)), where |a)i,2 is the 
generic element of 7ii,2- Using the inverse of the first of the (2.8), we obtain: 

MGi0) ly^ix) G-\e)\b)i,2 = MKi^)\b)i,2 . (2.21) 

Since the operator field z/g is defined on the Hilbert space Tig,^, eq.(2.21) shows that 
G~^{9)\a)i^2 is a vector of 7fe,/i, so G^^{9) maps Hi,2 to 7ie,n- G~^{6) : 7ii,2 ^ Ti-cn- 
In particular for the vacuum 1 0)1,2 we have (at finite volume V): 

\0),,,^G-\e)\0)^,2. (2.22) 

|0)e,/i is the vacuum for TYe,/*- In fact, from eqs.(2.8) we obtain the positive frequency 
operators, i.e. the annihilators, relative to the fields lydx) and ^'jn(x) as 

<:e^l,e = G-\9)ul^,al^,G{9), (2.23a) 

<l = G-\e) u\-^ al, G{e) , (2.236) 

vl7 Ke = G~\(^) /^M Gie) , (2.23c) 

PI, = G-\e) vl%* G{e) . (2.23rf) 

Eqs.(2.23) are obtained by using the linearity of operator G{6). It is a trivial matter 
to check that these operators do effectively annihilate \^)e,tM- 
Furthermore, for the vacuum state |0)e,/i the conditions hold: 

j (fx i/l(x)i/e(x)|0)e,^ = , j(fx vl{x)v^{xm,^^ = , (2.24) 



as can be verified using eqs.(2.8) and (2.22), or the definitions (2.23). 

In Section 3 we will explicitly compute eqs.(2.23), thus giving the dynamical map 
of the flavor operators in terms of the mass operators. 

We observe that G~^{9) — exp[9{S^ — S+)] is just the generator for generahzed 
coherent states of SU{2): the flavor vacuum state is therefore an SU{2) coherent 
state. Let us obtain the explicit expression for |0)e,^ and investigate the infinite 
volume limit of eq.(2.22). 

Using the Gaussian decomposition, G~^{9) can be written as [4] 

exp[0{S- — 5'+)] = exp{—tan9 S+) exp{—2ln cosO S3) expitanO S^) (2.25) 
where < ^ < f . Eq.(2.22) then becomes 

|0)e_^ = J\exp{-tane Sl)exp{-2ln cos9 S^) exp{tan9 5!)|0)i,2 . (2.26) 

k 

The right hand side of eq.(2.26) may be computed by using the relations 

53^0)1,2 = , 5i|0)i,2 7^0 , {Sim,,2^0 , (5i)3|0)i,2 = 0, (2.27) 

and other useful relations which are given in the Appendix A. The final expression 
for \0)e,n in terms of Sj_ and is: 



|0)e,M = n |0)e,M = n [1 + sin ^ cose {s1 -S'^) + l sin' 9 cos' 9 {{Str + {Slf 



+ 



sin' 9S'lS'l + -sin^9cos9(S^{S'l)' - S^S^)') + -sin" 9{S'l)'{S^)' 



|0)l,2 . 

(2.28) 

The state \0)e,n is normalized to 1 (see eq.(2.22)). Eq.(2.28) and eqs.(2.15) and 
(2.16) exhibit the rich coherent state structure of \0)e,tj,- 
Let us now compute i,2{0\0)e,n- We obtain 

i,2(oio)e,, = n (1 - ^ iMsxsmi,2 + \ Sin' 9 i,2msins'ir 10)1,2) (2.29) 

where (see Appendix B) 

iMsls'_\o),,2 =1,2 (o| EE«i«i,2)«2^:.,i)/5^.,i«^,2«i/5:Vilo)i>2 = 

a,T r,s 

rf „,s |2 



E I Ai<2 \'^Z,. (2.30) 



In a similar way we find 
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(2.31) 



Explicitly, Zfe is given by 



^ ^ /C^ [{uJk,2 + ^2) - (g^fc,! + mi)f 



where u!k,i = yk"^ + mf . The function depends on k only through its modulus 
and it is always in the interval [0, 1[. It has a maximum for k = ^Jn^ m2 which tends 
asymptotically to 1 when |m2 — mi\ — > 00; also, Zj. — > when k — > 00. 
In conclusion we have 

i,2(0|0)e,^ = n (1 - ^sin^^ Z,)' = nr(A;) = 

^-Qg«nr(fe) ^g^,znr(fe)_ (2.33) 

k 

Prom the properties of we have that T{k) < 1 for any value of k and of the 
parameters mi and 7712- By using the customary continuous limit relation Y,k ~^ 
/ d^k, in the infinite volume limit we obtain 

hm i,2(0|0)e,^= lim e(5^-^''''"'^^'^ = (2.34) 

y— >oo V— >oo 

Of course, this orthogonality disappears when ^ = and/or when mi = m2 (because 
in this case Zk = and no mixing occurs in Pontecorvo theory). 

Eq.(2.34) expresses the unitary inequivalence in the infinite volume limit of the 
flavor and the mass representations and shows the absolutely non-trivial nature of the 
mixing transformations (2.1). In other words, the mixing transformations induce a 
physically non-trivial structure in the flavor vacuum which indeed turns out to be an 
SU (2) generalized coherent state. In Section 4 we will see how such a vacuum struc- 
ture may lead to phenomenological consequences in the neutrino oscillations, which 
possibly may be experimentally tested. Prom eq.(2.34) we also see that eq.(2.22) is 
a purely formal expression which only holds at finite volume. 

We thus realize the limit of validity of the approximation usually adopted when 
the mass vacuum state (representation for definite mass operators) is identified with 
the vacuum for the fiavor operators. We point out that even at finite volume the 
vacua identification is actually an approximation since the fiavor vacuum is an SU (2) 
generalized coherent state. In such an approximation, the coherent state structure 
and many physical features are missed. 



3 The number operator, the dynamical map and the mass sectors 
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Wc now calculate the number of the particles condensed in the state |0)e,^. 

Let us consider, for example, the a^,! particles. As usual we define the number 

operator as iV^^ = J2r '^k\'^k i ^^'^ ^^^^ ^^^^ commutes with and with 

Sj. ioT k. Then we can write 



e,.(0|iV.MO)e,^ = e:,{0\KM,, (3-1) 

|0)g^ has been introduced in eq.(2.28). Prom eq.(3.1) and the relations given in 
Appendix C we obtain: 

e,M(0|^'jO)e,M = ^feSin2^. (3.2) 
The same result is obtained for the number operators N^^, Np^, Np^: 

e,^{0\N^p)e,^. = Zksin^e , a = a,P , i^l,2. (3.3) 

Eq.(3.2) gives the condensation density of the flavor vacuum state as a function 
of the mixing angle 9, of the masses mi and 1712, and of the momentum modulus k. 
This last feature is particularly interesting since, as we will see, the vacuum acts as a 
"momentum (or spectrum) analyzer" when time-evolution and flavor oscillations are 
considered. 

We remark that the eq.(3.2) (and (3.3)) clearly shows that the flavor vacuum |0)e,^ 
is not annihilated by the operators j, /5[ j, i = 1,2. This is in contrast with the 
usual treatment where the flavor vacuum |0)e/i is identifled with the mass vacuum 

|0)l,2. 

Notice that, due to the behaviour of for high k, expectation values of N^. 
are zero for high k (the same is true for any operator O^, i — 1,2, for which is 

l,2(0|Of |0)i,2 = 0) . 

In order to explicitly exhibit the dynamical map, eqs.(2.23), it is convenient to 
redeflne the operatorial parts of the flelds Uei^) and ly^ix) as (cf. eqs.(2.23)) w^' "o;^ ^ = 
%',e^fe,e> ^^^-i so ^^at write: 

^ G-\e) al, G{e) , (3.4a) 

al^^^G-\0)al,G{0), (3.46) 
f3l, = G-\d)Pl,Gid), (3.4c) 

f3l,,^G-\e)(3l,G{e). {3 Ad) 

We observe that and / = e,/x, depend on time through the time depen- 
dence of G{6). We obtain: 

^k,e = cosOal^^ + sin^ ^ («!i'«fc,2) + («M^^-fe,2) f^-k,'^ (3-5a) 
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^k,n = cos9al^ - sine 1] (KIs^i) + Kl2^^-fc,i) ^^-m) (^-S^) 

s 

= cose PL,,, + smeY: ((^1,2^:.,!) f^-k,2 + i<2vU,,) al],) (3.5c) 

s 

= cos^ PL,,, - smOY: ({v%,y_,,,) PL,,, + {ullvL,,,) all) (3-5^) 



Without loss of generality, we can choose the reference frame such that k = 
(0,0, |A;|). This implies that only the products of wave functions with r = s will 
survive (see Appendix B). Eqs.(3.5) then assume the simpler form: 



= cos0al, + sine ([/* a^, + V, P%) (3.6a) 



= cos e al,2 - sin 9 (u, al,, - V, PL^,,) (3.66) 

PL,,e = cose PL,,, + sine {U; PL,,, - V, al\,) (3.6c) 

PL,,^ = cos^ PL,,2 - sine {U, PL,,, + V, <J (3.6ci) 
with = (—1)'^ and 

U, = {ul\y,,,) = {v%,,vL,,2) (3.7a) 

V, ^ {ul],vL,,,) = -e^ iul],vL,,,) (3.76) 
where the time dependence of U, and V, has been omitted. We have: 

V, = \Vk\ e^^^^.^+'^fc.i)* , U, = \Uk\ e^('^'=-2-'^fe.i)* (3.8) 

\U,\^( ^^hl±^] ^ ( ->^,2 + mA ^ / ^\ 

\ 2uj,,, J \ 2u,,2 J \ {uJk,i + m,){cj,,2 + m2)J 

\V,\ = ( ^^J-l±^] " ( ^^hl±^] ' ( _ ^ (3.96) 
V 2u;,,, J \ 2u,,2 J \[uj,,2 + m2) {,ujk,, + m,)j 

\U,\' + \V,\^ = 1 , \V,\^ = ^Z, (3.10) 

For notational simplicity in the following we put Ui = uj,,i. 

It is also interesting to exhibit the explicit expression of |0)g^ in the reference 
frame for which k = (0, 0, |/c|) (see Appendix D): 

|0)e,M = n[(l-sin'^ \V,\^) - sine cose V,{A' + B'-)+ 

r 

+ e'' sin^ e V,{U;C'^ - U,D^) + sin^ e if A'-^'^l 1 0)1,2 (3.11) 
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with 



Al^al\(3% , Bl^al],P\, , Cl^a^n,! , Dl^al],n,2 (3-12) 

We observe that eqs.(3.6) can be obtained by a rotation and by a subsequent 
Bogohubov transformation. To see this it is convenient to put (cf. eq.(3.10)): 



— ^/t art 



and 



\Uk\ = cos©fe 



\Vk\ = sinGjfc 



so that eqs.(3.6) are rewritten as 



TT 



< e. < - 



(3.13) 
(3.14) 

(3.15a) 

(3.156) 

(3.15c) 
(3.15ci) 



where 



R^exphY: [Ki«^,2 + P%P-k,2) - K2«M + f^%f^-k,i) e"^^] (3.16) 



k,r 



k.r 



(3.17) 



^2 = exp\Yl ®fe 



fc,2 



(3.18) 



By use of these relations and noting that -R|0)i,2 = |0)i,2, we can "separate" 
the sectors {|O(0))i} and {|0(6))2} out of the full representation space {|0)e^} : 

{|O(0))i}0{|O(0))2}C{|O)e,,}. 

The states |O(0))i and |O(0))2 are respectively obtained as: 

|O(0))i = Sr'(0)|O)i = n H0fe + e'^' sin0fe /3:Vi«m) |0)i (3.19) 

k,r 

|O(0))2 = 52"'(0)|O)2 = n (cos0fe - e'- e^^^ sin0fc |0)2 (3.20) 

k.r 
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If one wants to work with the "mass" sectors {|0(G))i} and {|0(G))2}, the tensor 
product formahsm must be used, e.g. 

(Oi + 02) (|o(e))i <g) |o(e))2) = {Oi0i + i0O2) (|o(e))i |o(e))2) = 

= Oi|o(e))i ® |o(e))2 + |o(e))i ® O2|o(e))2 (3.21) 

with Oi, i = 1,2, any product of z/j neutrino field operators. For example, we 
have |0)i,2 = |0)i ® |0)2, al^^ = al^^ ® I, al^^ = I ® ^1,2, so that al^^al]^ = 
(«M ® ^) ® "m) = «fc,i ® <^k,2 and 

R-^al iR = cos e (a^^i (g) /) + e'''' sin ^ (/ (g) a]^ 2) (3-22) 

B2^R-^al^RB2 = cos ^ (a^^i ® /) + e'^ sin ^ (/ ® a^,2(0)) (3-23) 

with 

a^,,2(0) = cos e, al, + e'^^ sin 9, (3.24) 

and 

52"'|0)i,2 = (|0)i®|0(e))2) . (3.25) 

We note that |O(0))j, i = 1,2, are the vacuum states for Q;^j(0) — 
B{Q)-^al ^5(0), and (31^(0) = B{Q)-^(3l iB{Q), operators. 

We aiso note that , (|0)i ® |0(e))2) = , (|0)i (g) |0(e))2) = 0; 

«fe,^(|o(e))i®|o)2) = PlJ\m)i®m = 0, but V,^j|o(0))i® 10)2) ^ 0, 

|O(0))2) 7^ 0' etc.. Moreover, 2(O(0)|iV^f |O(0))2 = sin20fc, 
i(O(0)|iV,Y|O(e))i = sin'Qfc' and 2(O(0)|7V,Y|O(0))2 = i(O(0)|iV^Y|O(0))i = 0, 
a — a, (3, which show the condensate structure of the sectors {|O(0))i} , i = l,2. 
Finally, we observe that |O(0))i can be written as 

|O(0)), = exp (-^) \h) = exp (-^) \h) (3.26) 
with |/,) = exp (Efe,.(-l)^+'e^e^^'/3:V,«J) |0)i,2, and 

Sa, = - Y. ("fcli"fe,i log sin^ 0fc + "fe,i"fcli log cos^ 0fc) , i = 1, 2 . (3.27) 

k,r 

A similar expression holds for S|3^. It is known that (or 5*^.) can be interpreted 
as the entropy function associated to the vacuum condensate [7] . 

4 Neutrino oscillations 

We are now ready to study the flavor oscillations. In order to compare our result 
with the conventional one [1,2,5], we first reproduce the usual oscillation formula. 



12 



In the original Pontecorvo and collaborators treatment [5], the vacuum state for 
definite fiavor neutrinos is identified with the vacuum state for definite mass neu- 
trinos: \0)e,n = |0)i,2 = |0). As we have shown in the previous Section, such an 
identification is not possible in QFT; however, it is allowed at finite volume where 
no problem of unitary inequivalence arises in the choice of the Hilbert space. As 
already observed, even at finite volume, the vacua identification is only an approxi- 
mation. For shortness we refer to such an identification simply as to the finite volume 
approximation, meaning by that the approximation which is allowed at finite volume. 

The number operators relative to electronic and muonic neutrinos are 

= cos^ 9 al^^ial i + sin^ 9 a^^i'^\^2 + sin Q cos 9 {ol'^^^ol^^ + <^k,i'^\^\) , (4-1) 



= cos^ ^Q!fe|2'^fc,2 + sin^ ^Q^fe|i«fc,i - sin 9 cos 9{a'J^^al2 + «fc|2Q'fc,i) > (4-2) 
and, obviously, 

(0|Ar„':|0) = {0\N'^;\0) = 0. (4.3) 
The one electronic neutrino state is of the form 

K,e) = cos^ la^^i) +sin^ \ay. (4.4) 

where jctfe^e) = ctfelelO) ' l'^fe,i) = '^mI^) ; \'^k,2) = '^kl2\^)- The time evoluted of 
this state is controlled by the time evolution of \al i) and \cKl2) 

Keit)) = e-^"^^^'\al^) = cos^e--*|a^,,,)+sin^e-^-*|a^,,). (4.5) 
We have 

«el K: Ke) = 1 (4-6) 

and 

(«fc,e(^)l K'l Wl,e{t)) = cos^ ^ + sin^ ^ + 9 cos^ ^(g-^C'-^-'-i)* + e+*(^2-^i)i) _ 

= l-sin2 2^ sin^f— . (4.7) 



2 

The number of ctg particles therefore oscillates in time with a frequency given by the 
difference Aa; in the energies of the physical components ai and 0:2. This oscillation 
is a flavor oscillation since we have at the same time: 

KeWI WIM = sin2 2^ sin^ (^t) (4.8) 
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so that 

{aiM K: WIM + {^IM K: WIM = 1 • (4.9) 

Notice that the traditional derivation of eq.(4.7) is the same as the one presented 
above since (ci;^ g(t)| N^'J |ttfc,e(^)) — l(<^fe,el'^A;,e('^))P' "^^^ t)e seen from eqs.(4.4) 
and (4.5). Eqs.(4.7)-(4.9) are the well known results. 
Let us now go to the QFT framework. 

We have seen that |0)i,2 and |0)e,^ are orthogonal in the infinite volume limit. We 
choose to work in the physical (mass) representation |0)i^2 (same conclusions arc of 
course reached by working, with due changes, in the flavor representation \0)e,fi ) in 
order to follow the time evolution of the physical components |q;^_^) and |q;^ 2) ■ 

The number operators are now (cf. eqs.(3.6)) 

= cos' eallal, + sin' 9 pkf a^^a^^ + sin' 9 |14|' Pl2Pk,2+ 
+ sin^cos^ {u; all.al^ + Uk (x%al^^ + ^l^a^, + Vk all.Pl]^) + 

sin' 9 {Vk U, a^^^ + V,* U; dl^al,) , (4.10) 

= <A = 

= cos' 9all,al^, + sin' 9 |[/,|' a^al, + sin' 9 \Vk\' ^^,^1^+ 

- sin 9 cos 9 {U; a^al, + U, a^.^L - V; Pl^a^, - V, al],(3l\) + 

-e^ sin' 9 {V, f/* c^J^, + V,* U, Pl^a^,) . (4.11) 
The one electronic neutrino state is given by 

Ke) = <e|0)i,2 = cos^ + siu ^ C/fe [a^^) (4.12) 

with \al^) = a'i^\\0)i^2, i = 1,2. The action of the flavor operator a^^g is defined on 
|0)i_2 through the mapping (3.6) (see also eqs.(3.4)). Note that, as it must be, 

K,g|«y = 1 - sin'^ \Vk\' = 1 - 1,2(01 «^Xe I0)i,2 , (4.13) 

since ^ |0)i,2 7^ 0. 

The time evolution of the state |q;^ g) is given by 

Keit)) = e-^^.^*|a^,g) = cose e--*|ai;,i) + sin^ U, e^-^la^^^) • (4-14) 
We have: 

Ke\ Wle) = 1 - sin'e |y,|' = 1 - 1,2(01 Nl': |0)i,2 (4.15) 
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so that, 

Kel N'^: Ke) + 1,2(01 AT^^f |0)i,2 = 1 . (4.16) 

Wc note that the expectation value of N^^'^ in the vacuum 1 0)1,2 provides an essential 
contribution to the normalization equation (4.16) (see also eq.(4.13)). 

It is also interesting to observe that the term 1,2 (0| N^'J |0)i_2 plays the role of 
zero point contribution when considering the energy contribution of 0;^'^ particles. 

Note that we also have: 

i,2{0\Kf\0)i,2 ^ sm'e\VkW a^a,P, l^e,ii, (4.17) 

and 

«el Ke) = sin^^^ |\4f (1 - sin^^ \V,\') . (4.18) 

Eqs.(4.17) show the condensate structure of |0)i,2 in terms of definite flavor fields and 
of course are analogous to eqs.(3.3) which show the condensate structure of |0)e,^ in 
terms of definite mass fields. 
We also have 

= cos^^+|C/jt|' sin^^+|T4|' sin'^ cos" 9+\Uk\'' sm^ 9 cos" 9{e-'^'''-'^'^'+e+'^'''-'^'^') = 

- ( 1 - sin^^ \V,\^) - \U,\^ sm^29 sin^ (^^t^ . (4.19) 

This result reproduces the one obtained in the finite volume approximation (cf. 
eq.(4.7)) when \Uk\ 1 (and |Vjt| 0). The fraction of a^''' particles in the same 
state is 

KM^ZKM = 

= \Uk\^ sm^29 sin^ (^t^ + sin^^ l^fej^ (l - sin^^ , (4.20) 

where we recognize the contribution from the condensate in the state [a^ e) (cf- 
eq.(4.18)). Eq.(4.19) is to be compared with the approximated one (4.7). Note that 

KM WIM + KM WIM = Kel N'^: Ke) + K,e\ K,e) ■ 

(4.21) 

The normalization relation is written as 

KM K: KM+ 1.2(01 ivi^ |o)i,2 +KM ^Z KM - Kel Ke) = i ■ 

(4.22) 

We may also use Afk = {oik,e\<^k,e) ^ normalization factor and write eqs.(4.15), 
(4.18), (4.19), (4.20) and the normalization relation as 




(4.23) 
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^Kel^inO = sin^^lV^r (4.24) 
;^Ke(^)l ^al Keit)) = 1 - R, sm' 29 sin^ (^t) (4.25) 

^^Keit)\ ^'4 WkM = «in'2^ sin' (^^) + ;)^«el ^^a; l«^,e) (4-26) 

;^«e(^)l ^af Ke(^)) + ;^(4,e(^)l l<e(*)) " Ke) = 1- 

(4.27) 

respectively, with i?^ = ^ = i-liJem^ - 

In conclusion, eqs.(4.19) and (4.20) exhibit the corrections to the flavor oscillations 
coming from the condensate contributions. 

The conventional (approximate) results (4.7) and (4.8) are obtained when the 
condensate contributions are missing (in the |Vfe| — limit). 

Notice that the fraction of ctg''' particles which is condensed into the vacuum 
|0)i2 is "frozen", i.e. does not oscillate in time, as it is easily seen by noting that 
g-ii?i,2t|o)^ 2 = |0)i,2 (cf. also eq. (4.15)). 

It is remarkable that the corrections depend on the modulus k through |f/fcp = 
1 — \Zk. Since — >• when k ^ oo, these corrections disappear in the infinite 
momentum limit. However, for finite k, the oscillation amplitude is depressed by a 
factor \Uk\^: the depression factor ranges from 1 to | depending on k and on the 
masses values, according to the behaviour of the Zk function. It is an interesting 
question to ask if an experimental test may show such a dependence of the fiavor 
oscillation amplitude. 

We stress that the hmit A; — > oo of eqs.(4.19) and (4.20) gives an exact result and 
is not the result of the finite volume approximation. 

Since the correction factor is related to the vacuum condensate, we see that the 
vacuum acts as a "momentum (or spectrum) analyzer" for the oscillating neutrinos: 
neutrinos with k ^ y^mlrn^ oscillate more than neutrinos with k ~ ^/rnlfn2, due to 
the vacuum structure. Such a vacuum spectral analysis effect may sum up to other 
effects (such as MSW effect [8] in the matter) in depressing or enhancing neutrino 
oscillations. 

The above scheme is easily generalized to the oscillations in the matter. As 
well known [1,2,8], the two flavors oscillations picture is modifled due to the different 
interaction of and with the electrons of the medium. There are two contributions 
to this interaction: the first one, coming from neutral current, amounts to —GpTin/ \f2 
and it is present for both z/g and z/^; the second one, coming from charged current, 
is given by \f2Gpn(, and it is present only for z/g. Here and n„ are the electron 
and neutron densities, respectively. This produces a difference in the effective masses 
of i/g s-nd z/^, which can be expressed in terms of new free fields v\ and with new 
masses mi and -mi and a new mixing angle 9. 
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The mixing relations (2.1) are thus rewritten as: 

lydx) — i'lix) cos 9 + i>2{x) sin 9 

Uij,{x) — —ui{x) sm9 + D2{x) cos9 . (4.28) 

The tilde quantities are calculated by diagonalizating the hamiltonian of Ue and 
i/^ in matter, which is: 

jjnratter ^ ^~ _ Q^n^j^ ( 1 ) + ^^^""^ ( ) " ^^"^^^ 

The matter and the vacuum parameters are related as follows [1,2,8]: 

sin^2^~^sin^2^(^) (4.30) 
Atti^ = \{D - Am^ cos 29f + (Am^ sin 2^)^] ^ (4.31) 

\^(m\^m\ + D^ i^wf) (4.32) 

where D = 2^20 priek. 

When D = 2y/2GFnl"^k = Am^ cos 26' there is resonance and sin 26 goes to 
unity (MSW effect). 

In our scheme it is possible to treat oscillations in matter starting with mixing 
relations (4.28) and repeating all the procedure described above. So we can use all 
the results obtained in our treatment simply by substituting 9, mi, 1712 with the 
corresponding tilde quantities. In particular, the oscillation formula becomes in the 
matter: 

KM KM = ( 1 - ^^"^ \^k\') - \Uk\' sm'29 sm' (^t) . (4.33) 



and 



1,2 2 



5 Three flavors fermion mixing 

The extension of our discussion to three flavors is complicated by the proliferation 
of terms in the explicit computation of the quantities of interest. However, it is pos- 
sible to extract some results from the structure of the annihilators, without explicitly 
writing the expression for the vacuum state. 

Among the various possible parameterizations of the three fields mixing matrix, 
we choose to work with the following one: 

^ C12C13 S12C13 Sise"^ 

M= -S12C23 - Ci2S23Sl3e*'' C12C23 - Si2S23Sl3e*'' S23C13 | (5.1) 
V S12S23 - Ci2C23Sl3e"^ -C12S23 - Si2C23Si3e"^ C23C13 
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with Cij = cos 6ij, Sij = sin6'ij, since it is the famihar parameterization of CKM 
matrix [1]. 

To generate the M matrix, we define 

(^12 (6*12) = exp{9i2Li2) , (^23(^23) = exp{923 L23) , ^13 (^13) = exp{9i3Lis) 

(5.2) 

where 

L12 = J (^1(^)^2(3:^) — i^2(^)^i(^)) (5.3a) 

L23 = J - 4{x)t^2{x)) (5.36) 

Li3 = / d^x (e^' yl{x)u3{x) - e-'' vl{x)v^{x)) (5.3c) 

so that 

= G^2Gi}G23 G23GuGu (5.4a) 

ly'^ix) = G^,'GiiG,i u^ix) G23G13G12 (5.46) 

u^{x) = G^2'G^iG2i u^{x) G23G13G12 ■ (5.4c) 
The matrix M is indeed obtained by using the following relations: 

K(a;),Li2] = i^2"(a^) , K(a;),L23] = , [u'^ {x) , L,^] = e'' i^Si^) (5.5a) 
K(x),Li2] = -i^f(x) , K(x),L23]=i/3°(x) , K(x),Li3] = (5.56) 

K(x),Li2] = , K(x),L23] = -i^2"(^) , K(^),^i3] = -e-^'^r(^)- (5.5c) 

Notice that the phase S is unavoidable for three fields mixing, while it can be 
incorporated in the fields definition for two fields mixing. 
The vacuum in the fiavor representation is: 

\0)e,r = G^iG^iG2i\0)u3. (5.6) 

We do not give here the explicit form of this state, which is very complicated 
and is a combination of all possible couples al^iP^kj with i,j = 1,2,3. Nevertheless, 
we can obtain physical informations from the structure of the annihilators a^i, j3^i 
{I — e, fjL, t). In the reference frame k = (0, 0, |/c|) we obtain (see Appendix E): 

o^le = C12C13 + S12C13 {u^* + e'V^2 P-U) + e'' si3 (t/fa «^,3 + e'^13 P-is) , 

(5.7a) 

Qife,^ = (ci2C23 - e~'^ S12S23S13) Q;fc,2-(si2C23 + C""^ C12S23S13) (U^2 «M " ^''1^12 P-k,l) + 
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+ S23C13 {U^; al, + e^V^, , (5.76) 

"fe,r = C23C13 afc,3 - (ci2S23 + e~"^ S12C23S13) (t^23 "fc,2 " £^"^23 /^-U) + 

(5.7d) 

= (C12C23 - S12S23S13) /?:fe,2-(si2C23 + C12S23S13) (^^1^2 + 1^12 "m) + 

+ «23Ci3 (c/2'3 /3:.,3 - e'' ^2! «?3) , (5-7e) 

l3Lk,r = C23C13 /^Lfe.a - (C12S23 + S12C23S13) (C/23 P-k,2 + e'"1^23 ^fe^ + 

+ (^12^23 - e^' C12C23S13) (t/1'3 + e'^v,^3 ■ (5-r/) 

From eqs.(5.7) we observe that, in contrast with the case of two flavors mixing, 
the condensation densities are now different for different flavors (cf. eq.(4.17)): 

i23(0|7V,':|0)i23 = i23(0|iV^:|0)i23 = 1^121' + l^ill' , (5.8a) 
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123(0|iV,^;'|0)i23 = 

S12C23 + e~"^ C12S23S13 



123(0|iV«;|0):23 
2 



123(0|iV^;|0)i23 = 123(0|iV^nO)l23 



Cl2'S23 



'Sl2C23"Sl3 



|V^23 



k 1 2 
23 1 



'Sl2'S23 ~ e C12C23S13 



1^1 



k |2 
13 1 



(5.86) 



(5.8c) 



To study the three flavors neutrino oscillations, we observe that, in the finite 
volume approximation 



|afc,e) = C12C13 |afc,i) +S12C13 |afc,2) + e"^ Si3 |afc,3) 
The time evolution of this state is given by 



l«fc,e(^)) ^ C12C13 e "^'Vm) + *12Cl3 e' 

The number operator is 



-iu}2t I _,r 



fe,2/ 



+ e'" si3 e' 



-iuist I r 



fe,3/ 



(5.9) 



(5.10) 



rf r 



= C12C13 Qifc[iQifc,i + C12S12C13 a^JiQife^a + e'^ C12C13S13 oill^al^^+ 

+C12S12C13 Ofcls^M + Sl2C?3 «fcl2"fc,2 + S12C13S13 Q;fc|2"L3+ 

+e C12C13S13 afe;3"fe,i + e S12C13S13 afc;3"fc,2 + «i3 "fc'3"fe,3 • 



and thus 



Ke^ ^al Keit)) = 44 + 44 + 4 + 444(e 



„4 4 



4 „4 



i(a)2-wi)t _|_ g+i(a;2-a;i 



(5.11) 



2 2 2 /„-i(w3-W2)t I „+i(w3-W2)t\ 
•'12 ''13*131'^ ' ^ ) 



I.e. 



1 - cos^ ^13 sin^ 20i2 sin^ 



2 (UJ2— OJi 



t ] — sin^ 6*12 sin^ 26*13 sin^ 



2 a ■ 2csa • 2 / <^3 — ^1 
— cos 6^12 sm 2t/i3 sm I ^ 

On the other hand, the QFT computations give 

N'al Ke) = 1 " 44 1^2^ ' 4 1^1 
123(01 |0)i23 = 44 ln'2r + 4 I^1'3P , 



2 /UJ3 - L02 



k |2 
13 I 13 1 ' 



(5.12) 

(5.13) 
(5.14) 



and 
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+ {\U',,? inir + \U',,? \V^,? ) sl,clAz + \U',,? 444(e-'^'^^-'^^^* + e+'(---)0+ 

+ |C/^3|' 444(e-^('"3-'"i)*+e+^('^3-'^i)*)+|C/f2r It^fsl" 44s?3(e-^('"3-'"2)*+e+^('"3-'"^)*), 
namely 

-|C/f2r cos^^i3 sin2 2^i2 sin^ (^^^^^^) -|C/i'3l' 008^^12 sin^ 2^13 sin^ {p^t^ + 

-\Ui2? Wi^? sm^di2 sin2 2^i3 sin^ (^^^^l^t^ , (5.15) 
which is to be compared with eq.(5.12). 
6 Conclusions 

In this paper we have studied the fermion mixing transformations in the QFT 
framework. In particular we have considered the Pontecorvo mixing transformations 
for neutrino Dirac fields [5]. In the LSZ formalism of QFT [3,6,7] the Fock space 
of definite flavor states is shown to be unitarily inequivalent in the infinite volume 
limit to the Fock space of definite mass states. The fiavor states are obtained as 
condensate of massive neutrino-antineutrino pairs and exhibit the structure of SU (2) 
generalized coherent states [4]. The condensation density is computed as a function 
of the mixing angle, of the momentum modulus and of the neutrino masses. 

The neutrino oscillation formula is derived and its amplitude turns out to be 
momentum dependent. We suggest that such a result may be object of experimental 
investigation. 

In the k ^ 00 limit the momentum dependence disappears and the oscillation for- 
mula reproduces the usual one. Notice however that the oscillation formula we obtain 
in the limit A; ^ cxd is exact and is not the result of the finite volume approximation 
used in the conventional treatment. 

Since the oscillation term is related in our analysis to the vacuum condensate, the 
vacuum acts as a "momentum (or spectrum) analyzer" for the oscillating neutrinos. 
Such a vacuum spectral analysis effect may contribute in depressing or enhancing 
neutrino oscillations. 

Wc observe that the functional dependence of the oscillating term on the momen- 
tum is such that, if experimentally tested, may give indication on the neutrino masses 
since the function — 2 (1 — \Uk\'^) (cf. eqs.(3.10) and (4.19)) has a maximum at 
k — y/m^rn2. 
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Although the physically relevant quantities are given by expectation values of the 
observables, nevertheless it is interesting to observe that the ratio of the amplitudes 
of the |a;^ and [a^ 2) components of the state \cil^g{t)) is constant in time, as can 
be seen from eq.(4.14), and that such a feature persists even in the limit A; — > 00 (i.e. 
\Uk\ — > 1) where, however, the oscillation formula (4.19) reduces to the usual one 
(4.7). This is in contrast with the finite volume approximation where " decoherence" 
between the components and ja^ g) arises from the phase factor exp{—iAu! t) 

(see eq.(4.5)). 

We have shown that our discussion can be extended to the oscillations in matter 
and the corresponding oscillation formula is obtained. 

We have also studied the three fiavors mixing and have obtained the corresponding 
oscillation formula, which also in this case is momentum dependent. 

We stress the crucial role played by the existence in QFT of infinitely many 
unitarily inequivalent representations of the canonical anti-commutation rules. We 
have in fact explicitly shown that the neutrino mixing transformations map state 
spaces which are unitarily inequivalent in the infinite volume limit. In this way 
we realize the limit of validity of the identification of the vacuum state for definite 
mass neutrinos with the vacuum state for definite flavor neutrinos. We point out 
that such an identiflcation is actually an approximation since the flavor vacuum has 
the structure of an SU{2) generalized coherent state and it is only allowed at finite 
volume. The vacua identification is meaningless in QFT since the mass and the flavor 
space are unitarily inequivalent. 

Finally, we observe that although our discussion has been focused on the neutrino 
mixing, nevertheless it can be extended, with due changes, to other fermion mixing 
transformations, such as the CKM mixing transformations. In this last case, the so 
called free fields in the LSZ formalism are to be understood as the asymptotically 
free quark fields. 

The study of the QFT for mixing of boson fields with different masses is also 
in progress. In that case preliminary results show [9] that relations analogous to 
eqs.(3.2) and (3.6) hold so that we have a non trivial vacuum structure. In the boson 
case we find \Uk\ — coshcTfe and |Vjt| = sinhcTjt, with cr^ = ilogf^^) where Uki, 
i = 1, 2, is the boson energy. 
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Appendix A 

Using the algebra (2.19) and the relations (2.27), we have: 

^^^^io)i,2 = s'_si\o),,2 {sinsim,,2 = {smsim,,2 

{slfs'_\o),,2 = s1{siy\o),,2 + 2410)1,2 

{S^_rS'l\0),,2 = Sl{S'iy\0),,2 + 2S'_\0)^,2 

iSirSl\0),,2 = 6(4)2|0)i,2 (4)'5^|0)i,2 = 

{s^nsif \o)i,2 = 65^(5^)2|o)i,2 {s'_nsir\o),,2 = 6^^(5^)210)1,2 



{sins':m,,2 = 2A{sim,,2 
sis'_{sir\o)^,2^^sir\o),,2 

S|S'lSl\Q)^,2 = 
{S|rs1\0)^,2^{-lrs1\0)^,2 
(4)"4|0)i,2 = 5^10)1,2 

slsl{sim,,2 = S1{S1)'\0),,2 

Use of the above relations gives eq.(2.28). 

Appendix B 

Wave functions and Z^: 



{sins'im,,2 = o 
sisi{sir\o),,2^^sir\o),,2 
si{sms'im,,2 = o 

{S|r{S'_m^,2 = (-2)"(5^)2|0)i,2 
(4)"(4)^|0)i,2 = 2-(5^)2|0)i,2 
SlS>l{S'_f\0),,2 = -Sl{S'_f\0),,2 



u 



tr 



I 



a-k 



-r \ 



-,2t -2 



'^-k,l^k,2 — '^-k,l^k,2 



+ 



^a;jk,i + mi a;fe,2 + m2^ 
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-ki + iA;2 /ci — ik2 



Using the above relations eq.(2.32) is obtained. 
Eq.(2.31) follows if one observe that 

(^^)'|0)i,2 = 2 [{ul],v\,) {ull,v\,) - {ul],v\,) {ul],v\,)] all,P%^,aiy%\0)i,2 
Appendix C 

Relationships for the number operator and eq.(3.2): 

N^, Si] = E.,. ((Ol,2)<i«l,2 + K]ivuMli(^%) 

N^S'l] |0)i,2 = N^Sl\0)i,2 = N^{S'-)'\0)i,2 = 

N^,si\ si\o)i,2 = [n^,si] {siy\o)i,2 = [n^,si\ {simi,2 = 

N^, Si] |0)i,2 = 4|0)i,2 , [N^, Si] 5^|0)i,2 = (4)'|0)l,2 

N^Sl\0)i,2 = ^^|0)i,2 , 7V^(5^)2|0)i,2 = 2{Simi,2 



1,2 



1,2 



1,2 



1,2 



1,2 



1,2 



1,2 



1,2 



{0\S'iN^Sl\0)i,2 = Zk 



{0\S'lN^S'l{Simi,2 = i,2mS'iySlN^Sl\0)i,2 = {Zkf 
{Q\SXS^_NlSXS^_\Q)r,2^ Zk + \{ZkY 
{Q\{S^_fN^,{SXmi,2 = '2{Zkf 

k^^SXN^,St{Slf\Q)i,2 = Q{Zkf 
^SlNtSX{S^_f\Q)i,2 = 2 {Zkf 
^f{StfNlSXS>l_\Q)i,2 = iMSXS-N^{SXf{S^_f\Q)i,2 = 4 {Zkf 
msXf{StfNl{SXf{Stf\Q)i,2 = 24 . 



(0|(5^ 
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Eq.(3.2) is then obtained as 
e,,.(OK |0)e,^ = ^fc(sin2 e cos^ 6 + sin^ 9) + 

+ {Zkf sin^ 9{- cos2 ^ + 2 sin^ 9 cos^ 9-2 sin^ ^ + | cos^ ^ + i + § sin^ ^) = 
= Zk sin^ ^ . 



Appendix D 

For the calculation of |0)g^ it is useful to choose k = (0,0, \k\). In this reference 
frame the operators 5"^, 5"^, are written as follows: 



k,r 



k,r 

r 

k,r 

r 

k,r k,r 

where Uk, Vk have been defined in eqs.(3.7)-(3.10) and ~ {~^Y- It is easy to 
show that the su{2) algebra holds for S±^ and S^'^ , which means that the suk{2) 
algebra given in eqs.(2.19) splits into r disjoint SMfc,r(2) algebras. Using the Gaussian 
decomposition, |0)g_^ can be written as 

|0)e,^ ^X{exp{-tan9 S^/)exp{-2ln cos9 S^'"") exp{tan9 -S^''')|0)i,2 

r 

where < ^ < |. The final expression for |0)g^^ in terms of 5'^'' and S^''^ is then 
• M = n [l + sin ^ cos 9 [S^j'' - S'l!') - sin^ 9 S'l^^'S^j''] |0)i,2 , 



r 



from which wc finally obtain eq.(3.11). 
Appendix E 

Useful relations for three flavors mixing. 

We work in the frame k — (0, 0, and for simplicity we omit the k and the 
helicity indices. 
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GiiaiGi2 

G23 012G2Z 



Gi20i2Gi2 
G2Z (^^G22, 



, Gi2'^zGl2 



Gi2 0iG\2 
G 23 PIG 23 

GiiP2Gl3 
^^23^/33^^23 

GiiPsGis 
. GiiPlGu 



ai 



ci3 ai + e'^ si3 (c/i*3 as + V13 

C12 + S12 (U*2 012 + er V12 P2) 
C23 OL2 + S23 (t^2*3 0^3 + ^23 P3) 



G]_la2Gx3 = 0:2 



C12 01.2 
C23 Oi3 



Si2 (Uu ai 
S23 [Ui2 (y.2 



V,2 Pi) 
V,2 Pt) 



GiiasGis = ci3 as - e S13 (Uu ai - 1^13 



as 
Pi 



GTiPlGi3 = ci3 /3l + e'' s,3 {Ui3 PI - V^^ «3) 



= C12 /5l 

= C23 /32 



+ Si2 (c/12 - er V,*2 a2) 
+ S23 (U23 PI - ^2*3 «3) 



Pi 



C12 /3l 

C23 pI 
C13PI 

Pi 



Sl2 {U^2 Pi + er ^12 ttl 
•S23 



— e" 



Sl3 



with 



1 + 



{u;k,i + mi){u;k,j + mj)^ 
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1 1 



I '^i 



where i, j = 1, 2, 3 and j > i, and 

{U^2U^S - ^2*^23) = U^3 , (t^23^12 + t^l'* ^23) = ^ll 



k 
23 
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